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Abstract. We investigate soliton solutions of the Toda hierarchy on a quasi- 
periodic finite-gap background by means of the double commutation method 
and the inverse scattering transform. In particular, we compute the phase shift 
caused by a soliton on a quasi-periodic finite-gap background. Furthermore, 
we consider short range perturbations via scattering theory. We give a full 
description of the effect of the double commutation method on the scattering 
data and establish the inverse scattering transform in this setting. 



1. Introduction 

Solitons on a (quasi-)periodic background have a long tradition and are used 
to model localized excitements on a phonon, lattice, or magnetic field background 
(see, e.g., [S], [TT], [T3], [H], [T^, [H], [IB] and the references therein). Of course 
periodic solutions, as well as solitons travelling on a periodic background, are well 
understood. Nevertheless there are still several open questions. 

One of them is the stability of (quasi-)periodic solutions. For the constant so- 
lution it is a classical result, that a small initial perturbation asymptotically splits 
in a number of stable solitons. For a (quasi-)periodic background this cannot be 
the case. In fact, associated with every soliton there is a phase shift (which will 
be explicitly computed in Section 0]) and the phase shifts of all solitons will not 
add up to zero in general. Hence there must be something which makes up for this 
phase shift. Moreover, even if no solitons are present, the asymptotic limit is not 
the (quasi-)periodic background! A precise description of the asymptotic limit in 
terms of Abelian integrals on the underlying Riemann surface is given in [3] (see 
[lOj for a proof). In particular, the asymptotic limit can be split into parts, one 
which stems from the discrete spectrum (solitons) and one which stems from the 
continuous spectrum. 

The soliton part can be understood by adding/removing the solitons using a 
Darboux-type transformation, that is, commutation methods for the underlying 
Jacobi operators. Hence the purpose of the present paper is to complement [10] 
and provide a detailed description of the double commutation method when applied 
to a short range perturbation of a quasi-periodic finite-gap solution of the Toda 
lattice. In particular, we are interested in the effect of one double commutation 
step on the scattering data. 
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After introducing the Toda hierarchy in Section [2] and recalhng some necessary 
facts on algebro-geometric quasi-periodic finite-gap solutions in Section [3] we briefly 
review the single and double commutation methods in Section 3] and compute the 
phase shift (in the Jacobian variety) caused by inserting one eigenvalue for both 
methods. In Section [5] we review direct scattering theory for Jacobi operators 
with different (quasi-)periodic asymptotics in the same isospectral class. As our 
main result we give a complete description of the effect of the double commutation 
method on the scattering data. In addition, we provide some detailed asymptotic 
formulas for the Jost functions ip±{z,n) (which are normalized as n —^ ±oo) at the 
other side, that is, as rt — *■ =Fcxd. Our final Section[6]establishes the inverse scattering 
transform for this setting. Our main results here are the time dependence of both 
the scattering data and the kernel of the Gelfand-Levitan-Marchenko equation. 

2. The Toda hierarchy 

In this section we introduce the Toda hierarchy using the standard Lax formalism 
([!!])■ We first review some basic facts from [T] (see also [H]). 
We will only consider bounded solutions and hence require 

Hypothesis H.2.1. Suppose a{t), b{t) satisfy 

a{t) e 6(i) e £°°(Z,]R), a{n,t)^0, (n,t)eZxR, 

and let 1 1-^ (a(i), b{t)) be differentiable m e°°{Z) © £°°(Z). 

Associated with a{t), b{t) is a Jacobi operator 

(2.1) Hit):l^iZ)^e\Z), f^T{t)f, 
where 

(2.2) T{t)f(n) = a{n, t)f{n + 1) + a(n - 1, t)f{n - 1) + b(n, t)f{n) 

and £^ (Z) denotes the Hilbert space of square summable (complex- valued) sequences 
over Z. Moreover, choose constants cq = 1, c^, 1 < j < r, c^+i =0, set 

j 

9jin,t) = ^Cj-.i{5n,H{tf5n), 

(.=0 

j 

(2.3) hj{n,t) = 2a{n,t)'^Cj-i{Sn+i,H{tY6n) + Cj+i, 

e=o 

and consider the Lax operator 

r 

(2.4) P2r+2{t) = -H{ty+' + J2iMt)9j{t)S+ - hj{t))H{ty-^ + gr+i{t), 

3=0 

where S^f{n) — f{n± 1). Restricting to the two-dimensional nullspace Ker(r(t) — 
z), z e C, of T{t) — z, we have the following representation of P2r+2{t) 



(2.5) 



P2r+2{t) 



= 2a{t)GAz,t)S+ - Hr+i{z,t), 

Kcr(r(t)-z) 
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where Gr{z, n, t) and Hr+i{z, n, t) are nionic polynomials in z of the type 

r 

Gr{z,n,t) = ^z^gr-j{n,t), 

r 

(2.6) Hj.+i{z,n,t) = z''^^ + z^hj.-j{n, t) — gr+i{n, t). 

A straightforward computation shows that the Lax equation 

(2.7) ^i/(t) - [P2r+2{t), H{t)] =0, t e M, 
is equivalent to 

TLr{a{t),bit)), - d(t)-a(t)(5++i(i)-gr+i(0) =0, 

(2.8) TU(a(t), 6(t))2 - b{t) - {hr+i{t) - K+^{t)) - 0, 

where the dot denotes a derivative with respect to t and f^{n) = /(n± 1). Varying 
r e No yields the Toda hierarchy TLr{a,b) = (TL^(a, 6)i, TL^(a, 6)2) = 0. We will 
always consider r as a fixed, but arbitrary, value. 

We recall that the Lax equation (|2.7[) implies existence of a unitary propagator 
Ur(t,s) such that the family of operators H(t), t G R, are unitarily equivalent, 
H{t) = Ur{t, s)H{s)Ur{s,t). This also implies the basic existence and uniqueness 
theorem for the Toda hierarchy (see, e.g., [5D], [TH], or [31] Section 12.2]). 

Theorem 2.2. Suppose (ao,feo) € M = i°°{Z) © €°°(Z). Then there exists a 
unique integral curve t t-^ {a{t),b{t)) in C°°(K,M) of the Toda hierarchy, that is, 
TLr {a{t),b{t)) = 0, such that (a(0),6(0)) = (ao,5o). 

Finally, we recall the following result from 3J (compare also [20]), which says that 
solutions which are asymptotically close to a background solution at the initial time 
stay close for all time. 

Lemma 2.3. Suppose a{n,t), b{n,t) and a{n,t), b{n,t) are two arbitrary bounded 
solutions of the Toda hierarchy satisfying 112. 9\) for one to G M., then \2. 9]) holds for 
all t £ R, that is, 

(2.9) '^w{n)(^\a{n,t) -a{n,t)\ + \b{n,t) -b{n,t)\^ < 00, 
where w{n) > 0. 

3. Quasi-periodic finite-gap solutions 

As a preparation for our next section we first need to recall some facts on quasi- 
periodic finite-gap solutions (again see [1] or [21]). 

Let M be the Riemann surface associated with the following function 

23+1 

(3.1) <+2W = -nV^-^' Eo<Ei<---<E2g+u 

j=o 

where 5 € N and ^ is the standard root with branch cut along (— oo,0). M is a 

compact, hyperelliptic Riemann surface of genus g. A point on M is denoted by 

1/2 

p = (z, _^_2{z)) — (2, ±), z £ C, or p — oo±, and the projection onto C U {00} 
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by 7r(p) = z. The sets n± = {{z,±B}^g^^{z)) \ z e C\ Uf=o[-S2j, ^2j+i]} C M are 
called upper, lower sheet, respectively. 

Now pick g numbers (the Dirichlet eigenvalues) 

(3.2) (A,)f.i = (/i,,'^,)f=i 

whose projections lie in the spectral gaps, that is, fij € [E2j-i, E2j]. Associated 
with these numbers is the divisor I?^ which is one at the points /ij and zero else. 
Using this divisor we introduce 

zip, n, t) - [p) - (P^) - nA^_ {^+) + tU, ~ %^ e 

(3.3) z(n, f) = z(oo+, 71, t), 

where Sp^^ is the vector of Riemann constants, C/^ the ^-periods of the Abelian 
differential VLg defined below, and [OLpo) is Abel's map (for divisors). The 
hat indicates that we regard it as a (single-valued) map from M (the fundamental 
polygon associated with M) to C^. We recall that the function 0{z^{p,n,t)) has 
precisely g zeros flj{n,t) (with /ij(0,0) = fij), where 9{z_) is the Riemann theta 
function of M. 

Taking a stationary solution of TLg with constants Cj, 1 < j < g, as initial con- 
dition for another equation TL^. with constants Cj, 1 < j < s, in the Toda hierarchy 
(j2.8p one obtains the quasi-periodic finite gap solutions of the Toda hierarchy given 
by (see [211 Sections 13.1, 13.2]) 

2 ,2 %(n + l,t))g(z(n-l,0) 

"^^"''^ %(M)? ' 

d ( 9{w + z{n,t)) 



(3.4) 6^(„,t) = 6 + ^c,(g)— ln( 



_^ _ -jj \9{w + z{n—l,t))Jw=o 

The constants a, b, Cj{g) depend only on the Riemann surface (see [HI Section 9.2]). 
Introduce 

0(z{p,n+l,t)) ' rP 



(pq{p,n,t} ^ C{n,t) — — expl/ Woo+,oo_ 



(3.5) V',(p,«,i) = C(n,0,i)^J4^1^exp(n fcu^^^oo^+t f ils) 

e{z{p, 0,0)) V + J 

where C{n,t), C{n,0,t) are real-valued, 

,2 _ e{z{n - 1, t)) ^ _ 0{z{0, 0))e{z{-l, 0)) 



(3.6) C{n,t)'^ )-) ; , C{n,0,tY - 

u{z[n + l,t)} O{z[n,t))0[z_(n — l,t)) 

and the sign of C{n,t) is opposite to that of aq{n,t). Woo+,cxj„ is the Abelian 
differential of the third kind with poles at cx3+ respectively oo_ and Sis is an Abelian 
differential of the second kind with poles at cxd+ respectively cxd_ whose Laurent 
expansion is given by the coefficients [j + l)cs-j associated with TLs (see [2T| 
Sections 13.1, 13.2]). Then 

Tq{t)tljq{p,n,t) ^ 71(^)11) q{p,n, t), 

^'il}q{p,n,t) = 2aq{n,t)Gsip,n,t)ipq{p,n + l,t) - Hs+i{p,n,t)iljq{p,n,t) 

(3.7) = Pqas+2{t)^q{p,n,t), 
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where we use the hat to distinguish the quantities associated with TLs from those 



associated with TLg. 



The two branches ipq;±(.z,n,t) = tpq(p,n,t), p = (z,±), of the Baker-Akhiezer 
function are Unearly independent away from the branch points and their Wronskian 
is given by 



Here Wq{f, g) — aq{n){f{n)g{n+l) — f{n+l)g{n)) is the usual modified Wronskian. 

It is weh known that the spectrum of Hq{t) is time independent and consists of 
5 + 1 bands 

g 

(3.8) a{Hq)^[j[E2j.E2j+i]. 

For further information and proofs we refer to [21] Chapter 9]. 
Finally, let us renormalize the Baker-Akhiezer function 



(3.9) i^q{p,n,t) 



ijq{p,0,t) 

such that ipq{p, 0,t) — I and let us define as{p, t) via 

(3.10) exp(a,(p,i)) =V,(p,0,t) = C(0,0,t)^|=|^lMjexp(i 



4. Commutation methods and A^-soliton solutions 

In this section we investigate commutation methods when applied to a quasi- 
periodic finite-gap background solution. In particular, we compute the phase shift 
(in the Jacobian variety) introduced by the solitons. This can be found for the 
case of one-dimensional Schrodinger operators in [7] (see also [5] for the elliptic 
case) . The case of Jacobi operators seems to be missing and hence we provide the 
corresponding results to fill this gap. We want to be rather brief and refer to [S] or 
pn Ch. 11] for further details in this connection. Since the time t does not play a 
role in this section, we will just omit it. 

We start by inserting an eigenvalue using the single commutation method. Let 
Hq be a quasi-periodic finite-gap operator and let ipq^±{z,n) be the branches of the 
Baker-Akhiezer function which are square summable near ±oo. Fix Ai < inf a{Hq), 
cr G [—1,1], define 

(4.1) Uq^a{M,n) = i-t^?/;g,+ (Ai, n) + ^-^t/i, _(Ai, n), 

and let Hq a- be the (self-adjoint) commuted operator associated with 



(4.2) {Tq,„f){n) = aq,„{n)f{n + 1) + a,,,(n - l)/(n - 1) + bq,^{n)f{n), 
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where (see [211 Sect. 11.2]) 



Wg,o-(Al,n+ 1) 

M,,o-(Ai,n) Uq,^(Ai,7i + 1) 



bq^ain) = Ai - 05(71) 
(4.3) =6,(71) + 9* 



Ug,CT(Al,n+ 1) Mg^cr(Al,77) 
Q9('^)u9,a(Al,?l) 
U<;,(t(Ai, 77 + 1) 



iJq — Ai and Hqu — Ai restricted to the orthogonal complements of their corre- 
sponding one-dimensional null-spaces are unitarily equivalent and hence 



f^aciHq^a) — <^ac{Hq), CTsciHq^cr) — sc{Hq) 

{Ai}, (Te(-l,l) 
0, (7 €{-1,1} ■ 



Lemma 4.1. Let Hq be a quasi-periodic finite-gap operator associated with the 
Dirichlet divisor T^jx and letHq^a, —1 < cr <\, be the commuted operator associated 
with J/^.^l ). Then we have 

(4.4) aq.a{n) ^ aq,±i{n), bq^^{n) 6g,±i(7i) as n ^ ±00, 

where Hq,±i are the quasi-periodic finite-gap operators associated with the Dirichlet 
divisors T^fi^-^ defined via 

(4.5) «p„(I?^^^) = «p„(2?a) -^po(Pi) -^po(oo+)' Pi = (^i'±)- 

Proof. That i?q,±i is associated with the divisor T^[i^^ is shown in 211 Sect. 11.4] 
and the asymptotics follow since ^^ ^(Ai, n) ^ '3^Uq,±i(Ai, n) as 77 — > =poo. □ 

Similarly, one obtains the following result for the double commutation method. Let 
Ai e M\a^ss{Hq), define (see [211 Sect. 11.6, (2.30)]) 



1 " 



7 

j = ~oo 

(4.6) = i + iy,,„(7/.,,_(Ai),^,,_(Ai)) = -+V9,-(Ai,n)2</,,,_(Ai,77), 7^0, 
7 7 

and let -ff,,^ be the doubly commuted operator associated with 

^Cq^^[\i,n - l)cg,^(Ai,?i+ 1) 



a«,7("-) = a,q{n)- 
(4.7) 6,,^ (n) = 6,(71) -9' 



Cg,7(Al,n) 

a,(n)-0,^-(Ai,77)-0g^-(Ai,n+ 1) 



Cq,7(Al,n) 

Then 

Lemma 4.2. Let iJ, 6e a quasi-periodic finite-gap operator associated with the 
Dirichlet divisor Vj^ and let Hq^^, < 7 < co, be the doubly commuted operator 
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associated with i4-7\). Then we have 



aq(n)(l + 0(w(Ai)2") as n ^ -oo 

a,^ooN(l + 0(u;(Ai)-2») as n ^ +00 



(4.8) 



+ 0(w(Ai)-^") asn 



—00 



^ ~ \ 6,^ooH(l + 0(ii;(Ai)-2«) asn ^+00 ' 

where 10(2) = exp(Jp^' 1-^00+. oo_) is the quasi-momentum map and Hq^oc is the 
quasi-periodic finite-gap operator associated with the Dirichlet divisor defined 
via 

(4.9) apo(^A;,^) -fipo(^A) + 2Ap„(Ai), Ai = (Ai,+). 

Proof. Since the double commutation method can be obtained via two single com- 
mutation steps (see [HJ Sect. 11.5]), the result is a consequence of our previous 
lemma. The asymptotics follow from (14. 7p and p.Sp . □ 

Clearly, if we add k eigenvalues Ai, . . . , A^, then the asymptotics at +00 are given 
by the quasi periodic operator associated with the Dirichlet divisor 

k 

(4.10) «poP/i,, ;,J = «po(^Ai) + 2E^Po(^.)' A, -(A,,+). 

In particular, by choosing at least one eigenvalue in each gap, we can attain any 
prescribed asymptotics in the given isospectral class by Lemma 9.1 in [21j . 

Remark 4.3. If aq{n,t), bq{n,t) is a quasi-periodic solution of the Toda hierarchy 
and il)q{p,n,t) is the corresponding time dependent Baker-Akhiezer function, then 
aq^-y{n,t), bq,^{n,t) is a solution of the Toda hierarchy which is centered at 

(4.11) 2a(Ai)(n-v(Ai)t)+ln(7) = 0, 



wh' 



ere 



(Ai-) 1 /-(Ai.-) 



(4.12) a(Ai)=Re/ cjoo+,oo_, viXi) = r^Re ft,. 

5. Scattering theory 

In this section we review scattering theory for Jacobi operators with steplike 
quasi-periodic finite-gap background in the same isospectral class following [i] . Our 
only new result in this section will be a full description of the effect of the double 
commutation method on the scattering data in Lcmma l5.5l 

More precisely, we will take two quasi-periodic finite-gap operators associ- 
ated with the sequences a^, in the same isospectral class, 

a 

(5.1) aiH+) = a{H-) ^ S = |J [i?2„ i?2,+i], 

but with possibly different Dirichlet data {fi^}^j^i- We will add ± as a superscript 
to all data introduced in Section [3] to distinguish between the corresponding data 
of and . To avoid excessive sub/superscripts we abbreviate 

(5.2) Vq^(z,n) = Vq';±(^,"-) and Ti)q{z,n) = Tpq^^{z,n), 



8 I. EGOROVA, J. MICHOR, AND G. TESCHL 

that is, ^/j^{z,n) is the solution of decaying near ±00 and 'il!^{z,n) is the 
sohition of decaying near =Foo. Note that for A £ S we have 'ijj^{\,ri) — 
ii^t{\n). 

Let a{n), b[n) be sequences satisfying 

±00 

(5.3) 5]|n|(|a(n)-a±H| + |6(n)-6±(n)|) <oo 
and denote the corresponding operator by H . 

Theorem 5.1. Assume Ii5.3\) . Then there exist solutions 'il'±{z, .), z G C, of rip — 
zip satisfying 

(5.4) hin |«;(z)T"(^±(z,„)-V,^(z,„))| =0, 

where w{z) — exp{jl^^''^"' is the quasi-momentum map. 

Theorem 5.2. Assume i5.!^) . Then we have aess{H) = S, the point spectrum of H 
is finite and confined to the spectral gaps of , that is, <7p{H) = {pjjj^i C 
Furthermore, the essential spectrum of H is purely absolutely continuous. 

Using the fact that ip^{p,n) form an orthonormal basis for {dll+ , duj^) , where 

(5.5) duj^ = — -J- —dn, 

we can define 



±CXD 



(5.6) K±{n,m) = 2Re J ip±[\n)i}}^{\,m)du 
implying 

(5.7) ip±{z,n) = ^ K±{n,m)ip^{z,m). 

in—7i 

Next we define the coefficients of the scattering matrix via the scattering relations 

(5.8) i'^i^n) - a±(A)i/^±(A,n) +/5±(A)V'±(A,n), A e E, 
where 

,5..) o,,., = Z!fr'!t^ - IIfe#zf„,*_,A),*.(A„. 

hw = «^ ^ ^3k#^«'(.v(A),«Ay), 

1^(V'±(A),^±(A)) Rlil2{z) 

and Wn(f,g) — a{n){f{n)g{n + 1) — f{n + l)g{n)) denotes the Wronskian. Trans- 
mission T±(A) and reflection i?±(A) coefficients are then defined by 



o /3±(A) _ W^(V't(A),^±(A)) 



(5.10) T±(A) = a±^(A), i?±(A) 
The norming constants corresponding to pj £ ap{H) are given by 

(5.11) 7±i = El'^±(^J-'")l'' ^<J<1- 
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Note that 7±j- if pj coincides with a pole fi^ G n± of il:±[z,.). To avoid 
this, one could remove the poles by introducing V'±(^i-) ^ did in 4j. Since 
this normalization cancels out in the Gel'fand-Levitan-Marchenko equation and 
unnecessarily complicates the formulas below, we will allow zero norming constants. 
Moreover, ip±{Pj, •) = c^VtIPj' ■) '^i*^ ^ 

Lemma 5.3. The coefficients T±{X), R±{X) are bounded for A G S, continuous for 
A G S except at possibly the band edges Ej , and fulfill 



(5.12) T+(A)r-(A) + |i?±(A)|2 = l, AGE, 



(5.13) T±(A)i?±(A)+r±(A)i?^(A) =0, A G S. 
In particular, 

(5.14) |T±(A)p n ^:i4 + = 1' 

and hence |i?±(A)p < 1 with equality only possibly at the band edges {Ej}. The 
transmission coefficients T±{X) have a meromorphic continuation to C\S with sim- 
ple poles at fi^ if fi^ G Hip and simple poles at pj , 

R2g+2{P] 

nLi(Pj-Aif 

In addition, T±{z) G K. as z G M\E and 



(5.15) (Ites„r±(A))' = 7tj7-.J „, , _ 



^*(») = n ft 



""1 (") „=0 "9 (") 

The sets 

(5.17) S±{H) - {i?±(A),A G E; (p„7±j),l < J < g} 
are called left/right scattering data for H. 

Theorem 5.4. The kernel K±{n,m) of the transformation operator satisfies the 
Gel 'fand-Levitan-Marchenko equation 

(5.18) K±(n,m) + Y K±(n,l)F^{l,m) = ^(";'™) ±m > ±n, 

^ X±(n,n) 

w/iere 
(5.19) 

F±(m,n) = 2Re / i?±(A)^±(A, m)i^±(A, n)du;± + ^ 7±.,^,^(Pj , , m). 

TVote t/iaf i/ie apparent poles pf cancel with the zeros of duj^ and 'y±j at these 
points. 

The operator H can be uniquely reconstructed from S±{H) by solving the Gel'fand- 
Levitan-Marchenko equation. We refer to [4] for further details. 

Finally, we come to our principal new result in this section and investigate the 
connection with the double commutation method. The scattering data of the op- 
erators H, Hj are related as follows. 
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Lemma 5.5. Let H be a given Jacobi operator satisfying L5.3\} and choose Pq+i G 
M\a{H), 7 > 0. Then the doubly commuted operator H^, 7 > 0, defined via 
ip^lpq-^-i) as in Section^ satisfies 

(^20) a (n)'-^ [ "^^"^^ asn->-oo ^ / ^^(") as n ^ -oc 

[O.ZV) a^[n) I ^^^^^ as +00 ' \ b^{n) as n ^ +00 ' 

such that i5.3\) still holds, where is associated with Dirichlet data given by 

(5.21) apo(^Ai~) = £ipo(V) + ^AoiPi+^)- 

It has the scattering data 



%-(A,0))%+(A*,0)) ^ 



(5.22) 




.7(A) 


(5.23) 




.7(A) 


(5.24) 


T, 




(5.25) 


T_ 




where 






(5.26) 







0iz+{\,O)) 6l(z°°(A*,0)) 



B{\,pg+iyR+i\), 



l_ 6i(z°°(z,0)) 
C 9{z+{z,0)) 



B{z,Pq+i)T+{z), 
B{z,pq+^)T^{z), 



/ r{p-+) 
B{z, p) = exp - / 

\ Je(p) 

is the Blaschke factor and A = (A, +), z — (z, +). The constant C is given by 



..27^ / g(z°°((^+,0))g(z°°(cx3,,0)) 

^ ■ ' Y 0(z+(oo+,o))^^(z+(oo_,o)) ■ 

The norming constants corresponding to pj g ap{H), j — I, . . . , q, (cf. ^5.11]) ) 
remain unchanged except for an additional eigenvalue Pq+i with norming constant 
7_^g_l_i = 7. The norming constants 7+j^-y, j — I, . . . , q, are given by 

1 fe{z^{p,,o)r' 



(5.28) 7+j,7 = ^ (^ 6l(z+(p/o)) j 

anii t/ie additional norming constant ^+,q+i.y reads 
(5.29) 

^2 nf^i(p.+i -M;)Vgfe+(p;+i,o)) ^ , 

7+,9+i.7 = <^ ^ ^ or 00/-* ^T^+(pg+i)Res,=p,+,S(z,p,+i) 



7^29+2 (Pg+i) V^U°°(p;+i,0))" 
Proof. First we show that (|5.3p still holds. Note that 

a{n){l + 0{w{pq+,)^^) 



as n ^ — cxD 



""^^"^ ~ \ aoo(n)(l + 0(^«(p5+i)-2») as n ^ +00 

, , J 6(n)(l + 0(w(pq+i)2") asn^-cx) 
'^^"^ " 1 6oo(n)(l + 0(w(p5+i)-2") asn^+00 ■ 

Hence the asymptotics near —00 are clearly unchanged and for +00 it suffices to 
check 7 = 00. By Lemma 15.71 below . 

Coo{pq+l,n+ 1) ^ Coo{pq+l,n + l)ll!+{pq+l,n + 1) ^ C,^oo (Pg+1 , + 1) ^ ^^^^^ 
Coo(Pg+l,?^) Coo(p9+l,n)V'+(p9+l,n + 1) cJoo(Pg+l,n) 



SOLITON SOLUTIONS OF THE TODA HIERARCHY 



11 



where X^ngN ^ Therefore 

/r on^ I \ I ^ Vcoo(p9+l,n- l)Coo (Pg+1 , + 1) , > 

(5.30) aao{n)=a(n) r > a^,oo{n) 

such that X^riGN ~ '^q,oo{''T-)\ < and similarly for bao{n). 

Now we turn to the scattering data. By PU Lemma 11.19], the Jost solutions 
ip± n) of are up to a constant given by 
(5.31) 

C^(pg+l,n)Vj±(z,n) - — i— ■!/'-(Pg+l,n)W'„_l(V'-(Pg+l),'0±(2;)) 

u± .y(z, n) — ^ . 

Since this constant is equal to 1 for '0_^^(z,n) the fact that i?_ is unchanged 
follows from its definition and [U (11.107)]. The transmission coefficients are 
reconstructed from i?-(A) using 4, Theorem 3.6] and R+^j{X) follows then from 

i?+,,(A) = -^4^i?_,,(A). 

That the norming constants 7_j- are unchanged follows from f2T, Lemma 11.14]. 
For 7+j-,7, j = 1, . . . , g + 1, we use (15.15^ 

_nLi(P.-M,^)'(ReSp^.T±,,(A))2 



7+,. 



i?2s+2(Pi)7-J 

□ 



Remark 5.6. // we choose p — pj £ ap{H) and 7 — — 7_ j, then the eigenvalue 
Pj is removed from the spectrum and it is straightforward to see that an analogous 
result holds. 

The following result used in the previous proof is of independent interest. 

Lemma 5.7. Let H be a given Jacobi operator satisfying I15.3\) . Then for every 
z G C\S] and every k G N we have 

(5.32) i}^{z,n)i}±{z,n + k)-a±{z)tp^{z,n + k)'4}^{z,n) = {k+l)\w{z)\''C±{z,n), 

where I]±„gN "-|C'±(z, ")| < 00. 
Similarly, one has 

(5.33) Cooiz,n)ip±{z,n + k) - a±{z)c^^^{z,n)%l)'^{z,n + k) = C±{z,n,k), 



where X]±neN "'l^±(-^' ^)l < °^ '^"'^ Coo{z,n), c^^^(z,n) are defined as in (4- 

Proof. The proof is an extension of [14[ Lemma 3.4] and we will only consider the 
'+' case. First recall that the Green's functions of iJ+ and H are given by 

X il^+iz,m)^+{z,n) ij^{z,m)^+{z,n) 

ujlz.n.m) — — - — — i , G[z,n,m) — n > m, 

' ^ M^,(^+(z),V'|(z)) ' - 

respectively. Considering matrix elements in the second resolvent identity (H — 
z)-^ - {H+ - z)-i ^{H - z)-^{H+ - H){H+ - z)-i we obtain 

C+{z,n) = ^^jj^ E W{iP-{z),i,+ {z))G{z,n,m){H+ - H)G+iz,m,n + k) 
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for z e C\a{H). Since the poles of W {■t/j-{z) , ^p+{z))G{z , n, m) at z e (Jp{H) are 
removable, the formula holds for all z E C\I]. Estimating the right hand side using 
\G{z,n,m)\ < const I and \G+{z,n + k,m)\ < const we 

obtain 

|C+(z,n)| < ^^^^ E kWll—l+l-^+'llalaM-a+MI + lbM-b+MI). 

We split the sum into three parts 

( |w(z)|'=^|w(z)|2|"-™l m < n, 
|^(^)||«-m|+|„-m+fe| ^ ) \w{z)\'' n<m<n + k, 

[ |w(z)|'=|u;(z)|2|"-'"+'=l m>n + /fc. 

Since (15. 3p holds for c(m) = 2|a(rn)— a+(m)| + |6(m)— 6+(m)| and we have |w(2:)| < 1 
for z e C\E, we can apply Lemma lA. II to verify X^nGN < 
The claim (|5.33p is a consequence of (|4.6I) and (|5.32p 

n 

j = -oo 

a+(z)^c+^(z, n)tp+{z, n + /c)^ + C'±(z, n, k)%p'^{z, n + A:), 
where C±(z, n, /c) again can be estimated using Lemma lA. II □ 

6. Inverse scattering transform 
Let a{n,t), h{n,t) be a solution of the Toda hierarchy satisfying 

±CXD 

(6.1) E|n|(|a(n,t)-a±(n,t)| + |fe(n,t)-6±(n,t)|) < (X3. 

n=0 

Note that by Lemma [231 it suffices to check (|6.ip for one to G IK (as background 
take H~{t) and insert g eigenvalues such that the asymptotics on the other side 
are given by H^{t)). 

Jost solutions, transmission and reflection coefficients depend now on an addi- 
tional parameter t e R. The Jost solutions '>p±{z, n, t) are normalized such that 

V'±(z, n, t) = 'ipfiz, n, t) (1 + o(l)) as n ±00, 



where (c.f. ([SJ]) ) 

(6.2) V',^(z,n,t) = ' ' =:exp(-af(z,t))^±(z,n,t). 



Moreover, we set 

(6.3) exp(aj(z,t)) = -0^(2,0, t). 



Note that we have exp(Q;^ (z, t)) = exp(a^(z,t)) for A € S. 

Transmission and reflection coefficients are then defined via the normalized Jost 
solutions V-'± (2, n, t). Moreover, 

(6.4) <y{H{t)) = (j{H), t G M. 

To avoid the poles of the Baker- Akhiezer function, we will assume that none of 
the eigenvalues pj coincides with a Dirichlet eigenvalue /.t^(0,0). This can be done 
without loss of generality by shifting the initial time to = if necessary. 
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Remark 6.1. Due to this assumption there is no need to remove these poles for the 
definition ofj±j, as we did m [2], [I]- Since the Dirichlet eigenvalues rotate in their 
gap, the factor needed to remove the poles would only unnecessarily complicate the 
time evolution of the norming constants. Moreover, these factors would eventually 
cancel in the GeVfand-Levitan-Marchenko equation, which is the only interesting 
object from the inverse spectral point of view in the first place. 

Lemma 6.2. Let {a{t),b{t)) be a solution of the Toda hierarchy such that k5.3^) 
holds. The functions 

(6.5) %l'± [z, n, t) = exp(Q;^ (z, t))rp± (z, n, t) 
satisfy 

(6.6) H(t)ip±{z,n,t) ^ zip±{z,n,t), ^ip±{z,n,t) = p2s+2{t)il^±{z,n,t). 

Proof. We proceed as in [3], [20l Theorem 3.2]. The Jest solutions tp±{z,n,t) 
are continuously diflFerentiable with respect to t by the same arguments as for z 
(compare [U Theorem 4.2]) and the derivatives are equal to the derivatives of the 
Baker-Akhiezer functions as n — > ±cxd. 

For z G p{H), the solution u±{z, n, t) of (|6.6p with initial condition ip±{z, n, 0) S 
£j_(Z) remains square summable near ±oo for all t S M (see [19] or [211 Lemma 
12.16]), that is, u±{z, n, t) = C±{t)ip±{z, n, t). Letting n ±oo we see C±{t) — 1. 
The general result for all z G C now follows from continuity. □ 

This implies 

Theorem 6.3. Let {a{t),b{t)) be a solution of the Toda hierarchy such that (5. 
holds. The time evolution for the scattering data is given by 

(6.7) T±{z,t) = r±(z,0)exp(aj(z,t) - a^{z,t)), 

(6.8) i?±(A,<) = i?±(A,0)cxp(a^(A,t) -a^(A,i)), A G E, 

(6.9) 7±j(t) = 7±,,(0)exp(2a±(p,,t)), I < ] < q. 



Proof. Since the Wronskian of two solutions satisfying (|6.6[) does not depend on n 
or t (see [19], [21] Lemma 12.15]), we have 



^ W{4,±{z,t),i,^{z,t)) ^ exp(aT(^^t) + a,±(z,t)) W{4,±{z,t),i,±{z,t)) 

W{ij^{z,t),ij±{z,t)) e-x:p{at{z,t) + af{z,t))W(^^{z,t),i}±{z,t)) 
= exp(aj(z,i) - af{z,t))T±{z,Q). 

The result for R±{X,t) follows similarly. The time dependence of j±j{t) follows 
from U±{pj,.,t)\\ = \\Us{t,0),p±ip,,.,0)\\ = l|^±(p,, .,0)1|. □ 



Corollary 6.4. The quantity T±{X,t)T^{\,t) = 1 - \R±{X,t)\'^, A G E, does not 
depend on t. 

Another straightforward consequence is: 
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Theorem 6.5. The time dependence of the kernel of the Gel'fand-Levitan-Marchenko 
equation is given by 

F^{m,n, t) ^2Re J i?±(A, 0)V'^(A, m, t)V'^(A, n, t)dw±(0) 

+ X] T±.J (O)V'^ (Pj > V"^ {Pj ,n,t). 
i=i 

Proof. Just employ Theorem 16.31 to rewrite 

F±(m,n,t)= / R±{p,t)^f{p,m,t)^^{p,n,t)duj^{t) 

+ (^)'^q^ ^n,t)4}^ {pj ,m,t), 

where we also use that exp(a^(A, t) + af{X, t)) = Ilj^i a-ai^(o'o) ' 
Finally we note ([IS], [H Section 14.5]) 

Lemma 6.6. Let {a{t),b{t)) be a solution of the Toda hierarchy such that i5.3\) 
holds. Choose p G M\(t{H) and 7 > 0. Then (a-y{t),b-y{t)) defined via ip-{p,n,t) 
using the double commutation method is again a solution of the Toda hierarchy such 
that i5.S\) . with H^(t) accordingly changed, holds. 

Acknowledgments. I.E. gratefully acknowledges the extraordinary hospitality 
of the Faculty of Mathematics of the University of Vienna during two stays in 2006, 
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Appendix A. Some estimates 



In the proof of Lemma 15.71 we need the following elementary result. Consider a 
sequence c € and abbreviate 

00 00 
(Al) ||c||^=sup|c(n)|, ||c|li = ^|c(n)|, ||c||ia = ^ n|c(n)|. 

n— n—1 

Lemma A.l. Suppose w is some complex number with \w\ < 1 and c € £(7j) 
satisfies ||c||oo, ||c||i,i < 00. 
Then 



and 



00 

||^c(n + m)u;"||i.i<^-^||c||i,i 



E°° / N mil 1 II II M II II IwP 

c 71 - mK" 1,1 < 1 ^ c 1,1 + ■7T^T^Mi + 



SOLITON SOLUTIONS OF THE TODA HIERARCHY 



15 



Proof. The first estimate follows from 

oo oo oo 

II ^ c(n + m)w;'"||i,i = X! "| ^'^('^ + 'm)w"' 

m=0 n—1 m— 

oo oo 

- ^^{n + m)\c{n + m)\\wr 

n=0 m=0 

oo ^ 

m=0 ' ' 

Similarly, the second follows from 

oo oo oo 

II E c(n - m)w"^\\i^i = E "I E c(n - m)w"^ 

m=0 n— m— 

E°° / m(m — 1) 
I 2 11^11°° + "^l|c||i + 

m=0 ^ 

□ 
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